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We present a time-domain analysis of carrier dynamics in a 
semiconductor superlattice with two minibands. Integration 
of the density-matrix equations of motion reveals a number 
of new features: (i) for certain values of the applied static 
electric field strong interband transitions occur; (ii) in static 
fields the complex time-dependence of the density-matrix dis- 
plays a sequence of stable plateaus in the low field regime, and 
(iii) for applied fields with a periodic time-dependence the dy- 
namic response can be understood in terms of the quasienergy 
spectra. 
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Bloch oscillations (BO) are one of the most striking 
predictions of the semiclassical theory of electronic trans- 
port: in any system of independent electrons in a periodic 
potential the electron velocity becomes a periodic func- 
tion of time with characteristic frequency ujb = eEd/Ti, 
where d is the lattice period and E is the applied field 
. In ordinary bulk materials these oscillations cannot 
be seen, because collisions dephase the coherent motion 
of electrons on a time-scale which is much shorter than 
Tb = 2'k/u)b- However, as pointed out by Esaki and Tsu 
the conditions for observing BO's are much less strin- 
gent for high-quality semiconductor superlattices. Re- 
cent years have witnessed an intense experimental activ- 
ity in this area, culminating in the observation of tera- 
hertz radiation from coherently oscillating electrons Q]. 

There has been equal activity on the theoretical side. 
Holthaus 1^ analyzed the semiclassical motion of elec- 
trons in a single mini-band subjected to a strong alter- 
nating electric field. Studies of this kind have gained im- 
portance due to the emerging free-electron lasers, which 
open the possibility of experimental probing of the theo- 
retical predictions. For certain values of the system pa- 
rameters a dynamical localization takes place Q: the 
average velocity vanishes. This phenomenon can alter- 
natively be called band collapse Ignatov et al. 0] 
pointed out an interesting analog between Bloch elec- 
trons and the Josephson effect. Very recently, Meier et 
al. considered coherent motion of photoexcited car- 
riers in the presence of Coulomb interaction, and found 
out that BO's should persist even in the limit where the 
exciton binding energy is comparable to the miniband 
width. 

The papers quoted above have mainly concentrated 
on studying systems with one miniband the central 



theme in the present work is to study the dynamics of 
electrons in a two-hand superlattice. The second mini- 
band adds an essential feature to the model: it is possible 
to study how Zener tunneling affects the dynamics of the 
carriers. Our method consists of setting up, and solving, 
the density-matrix equations of motion for the two-band 
system. In this paper we focus on the coherent part of the 
motion. This coherent motion displays in its own right 
a number of interesting features, which we shall describe 
after having sketched the general formalism. 

First we need to define the microscopic model under- 
lying the density matrix calculation. The model Hamil- 
tonian is 

H = Y^ [(^S + neEd)aian + (A^ + neEd)blhn 

n 

— ^(4+ia« + 4a„+i) + -^(p^n+iK + hlb^+i) 
+eER{alX + bla,,)] . (1) 

The integers n label the lattice sites and the operators a 
and h refer to electrons in the two minibands; the first 
two terms give the (field-dependent) site-energies, the 
next two describe site-to-site hopping, and the last one 
is the term responsible for the interband transfer. The 
overlap matrix-element R is model-dependent and may 
depend on time through its momentum argument; we 
take it as a constant corresponding to a Kronig-Penney 
model. At zero applied field (|l|) leads to two mini- 
bands, e°^''(fc) = Ao''' T (Ai'V2)cos(A:d), while at finite 
static fields the spectrum consists of two interpenetrat- 



ing Stark-ladders |10| 



Ag- 



' eEdr" " - meEd. The 



numbers r"" [E) must be determined numerically, and 
results of such a calculation are shown in Fig. ^ It is 
important to notice that for certain field values the two 
levels come very close to each other; as we shall show 
below this leads to profound effects on Zener-tunneling. 

Let us next consider the equation of motion for the 
density matrix for a- and 6-electrons. With accelerated 
Bloch states as the basis [|ll|, the diagonal elements of 
the density matrix give the electron density at a /c-point 
following the semiclassical trajectory in reciprocal space, 
k{t) = K ~ [e/h) y^Eit')dt'. We assume equilibrium 
at t = 0, when the density matrix is diagonal. Defining 
P± {K, t) = p%- (t) ± {t) , we find the following equations 
of motion lfl2| : 



p+{K,t)=0 



(2) 



1 



-{K,t)= -Re\h{K,t) dt'h*{K,t')p-{K,t')^ 



(3) 



where 



h(K,t) = 2-E(t)Rcxp 

h 



/ Ae{k{t'))dt' 



(4) 



In the above equations we have defined Ae(fc) = — 
e°(fc), and assumed that the intraband coupHngs are 
identical |13[ . 

Eqs. (0-0) require several comments. If the interband 
coupling is turned off, they reduce to the normal col- 
lisionless Boltzmann equation (for two minibands), and 
thus contain, as special cases, the following standard re- 
sults: (i) for static field one finds standard Bloch oscilla- 
tions; and (ii) a harmonic £'-field leads to the analog with 
Josephson effect of Ref. Q], and in particular, reproduces 
the band collapse discussed by Holthaus Note that 
only p_ is affected by interband transitions, while p+, 
which fixes the particle density, is a constant of motion. 

By differentiating the equation of motion for with 
respect to time, one finds [O 



u . 

-P- 



+ u^p- +u(t) I dt' 



= . 



(5) 



Eq.(||) forms the basis of our analysis, and the rest of 
this paper will describe the numerical results obtained 
from it under a number of specific physical conditions. 

a. Steady fields. In this case the dc-field is turned on 
abruptly at < = 0. From Eq.(^) u is time- independent, 
and Eq.(|5|) can be reduced to an ordinary third order 
differential equation: 



= 



(6) 



The accompanying initial conditions are P-{0) = p'L, 
P-{0) = 0, and p_(0) = -uV-(O). Thus, the ini- 
tial values are determined by the miniband parameters, 
the fc-point in the Brillouin zone and the temperature 
(which enters through p_{0)). From Eq.(||) it follows 
that p^{K,0) can be chosen equal to 1 without loss of 
generality. 

Fig. H displays a typical time-dependence of p_, ob- 
tained from Eq.(^) by numerical integration. One ob- 
serves a very sensitive behavior with respect to variations 
of the applied field: for certain field values an 'inversion' 
takes place: p_ reaches —1, which is the negative of its 
initial value, while for other nearby field values p_ stays 
close to its initial value. This behavior can be under- 
stood by examining Fig. ^ the two energy levels are very 
close to each other for certain electric field strengths, and 
in the corresponding neighborhoods a strongly enhanced 
band-to-band transfer takes place. This situation is quite 
different from what one expects from simple Zener tun- 
neling theory j^] : there the tunneling rate is a monotonic 



function of the applied field. The situation is summa- 
rized in Fig. ^, where we plot these Zener resonances as 
a function of the applied field. One can assign an in- 
dex to the resonances: the resonance at the highest field 
(which corresponds approximately to aligning the levels 
at adjacent quantum wells), is called the first resonance, 
the next highest the second resonance (the case of Fig. 

and so forth. Adopting this numbering scheme we 
observe that in the low field limit i?'"^ « A^°'/ned (here 
A'"^ = - Ag). 

Further insight into the physical meaning of the vari- 
ous oscillations of Fig. ^ can be obtained by considering 
the Fourier transform of Let us first try to estab- 
lish a qualitative picture of what to expect. The initial 
state of the system is described by some wave-function, 
say |\E'(0)), which is not an eigenstate after the field has 
been turned on. However, it can be expanded in terms 
of the eigenstates: |5'(0)) = X)™ "^"^n ~^ dni^n- Since 
p{t) = |^(t))(4'(i)|, and each eigenstate evolves accord- 
ing to i}n ''{t) = exp(— ie^'^'i)?/'^''', large Fourier compo- 
nents in p_ are expected to occur at hto — meEd and 



Hlu ^ ±eEd{r°- - r'') + meEd |14|. This expectation 



is fully born out by the numerical evaluation of p-{uj). 
Fig. U shows the results of the two independent calcu- 
lations: the continuous lines are obtained based on Fig. 
|l|, while the asterisks come from the Fourier transform 
of p- (t) . Naturally, the more laborious calculation based 
on p-{t) contains also more information: the magnitudes 
of the Fourier components are needed in the evaluation 
of other physical quantities, such as the current, which 
will be addressed elsewhere | p^ . 

It is also of interest to examine the effect of vary- 
ing the superlattice parameters. Fig. ^ shows the time- 
dependence of p-{K = 0) when the field is tuned to 
the eighth Zener resonance, and we have increased the 
bandwidths and interband coupling. A distinctive set of 
stable plateaus has developed. The transitions between 
the plateaus occur at instants t — ■^Tb, fTs, fTs--- after 
the field was turned on. Thus, the life-time of a plateau 
is (approximately) equal to the Bloch period, and transi- 
tions occur every time the fc-point reaches the Brillouin 
zone edge ||l5[ . This behavior is generic to the low field 
regime, E < A^^/ed, and we can qualitatively under- 
stand features in the time-dependence of p_ {t) by consid- 
ering the semiclassical motion of a fc-point between the 
extrema of the Brillouin zone; transitions to the other 
miniband occur mainly at zone-edges, where the energy 
separation between the minibands is at minimum. 

We can also understand the number of oscillations on 
a given plateau by examining Eq.(|^) under some sim- 
plifying assumptions. In particular, if we assume that 
Ae has a weak time-dependence, it is easy to solve (|^) 
analytically. The solution suggests defining a 'local' 
time-dependent frequency for a general, but sufficiently 



slowly varying Ae to be (i) 
LOc = 2(|i?|/c?)a;B and hujdik) 



Lo^ + Luj{k{t)). Here 
Ae(fc). Thus, in the 



2 



low field limit we can identify the number of oscillations 
A^osc = {i^i{t))/uiB = A^'^/Uujb, where the time- average 
was calculated over the Bloch period Tb- Consequently, 
at the n:th Zener resonance, we find A'osc = n . In 
Fig. H one can distinguish two periods of oscillation in 
any of the plateaus (even though the plateaus are not 
very clearly resolved for this particular set of parame- 
ters) , while Fig. || clearly shows eight periods of oscilla- 
tion within a plateau. 

In the high field regime, E > A^"/ed, the situation 
differs drastically from semiclassical expectations: the 
plateaus vanish, and we find from Eq.(^), both numer- 
ically and analytically, that p^{t) oscillates between -1 
and -t-1 with a single frequency ojc- 

h. Alternating fields. In the case of a temporally 
periodic driving field one can make a close parallel to 
the treatment of a spatially periodic potential. Given a 
Hamiltonian with H{t + T) = H{t), the starting point is 
Floquet's theorem applied to the Schrodinger equation. 
This gives wave-functions of the form 



ipc = cxp(— iei)ue(i) , 



(7) 



where u^it + T) = u^{t). We have introduced the 
quasienergy e, which is defined modulo ujac = 27r/T, lead- 
ing in a natural way to the definition of a quasienergy 
Brillouin zone. The functions are eigenstates to 
S = H ~ ih-^, i.e. they are given by the temporal 
Bloch theorem applied to S. These considerations have 
been employed by Holthaus Q in his analysis of the one- 
miniband case, and we now wish to extend these con- 
cepts to two minibands in the context of the Hamiltonian 
(|I|) . The quasienergy spectrum for an applied field of the 
form E{t) — El cosujact is shown in Fig. ^ For van- 
ishing band coupling the band collapses occur whenever 
eEid/hiVac equals a root of the zeroth Bessel function 
(the corresponding field strenghts are indicated by ver- 
tical lines in Fig. |6|), and we observe that the two-band 
model displays strict band collapse at only one of these 
roots. However, there is a clear tendency towards band- 
width narrowing at the other roots. Thus, band collapse 
seems to be a generic feature, and not a special feature 
of the previously studied one-band models. 

The next step in the analysis is to solve the equation 
of motion for /9_ for the time-dependent case; one must 
now use (|^) as a starting point. One can perform a com- 
pletely analogous construction as is done in Fig. ^; now 
the quasi-energy differences (modulo iVac) correspond to 
the leading Fourier components of p_ (w) . There is an im- 
portant difference between the static and time-dependent 
cases : in the static case the energy spectrum does not 
depend on how the field reaches the final value, while in 
the time-dependent case the quasienergy bands are trans- 
lated in reciprocal space by an amount which is the neg- 
ative of the semiclassical change in crystal momentum 
during the turn-on period. The existence of the corre- 
sponding frequencies in the Fourier spectra of p- can also 



be directly derived from the equation of motion for p- 
[|l2|. In the low field fimit, Ei < A'"'/ed, we find that in- 
terband tunneling is quenched, and the time-dependence 
is described by the semiclassical equations of motion. 

In summary, we have presented a time-dependent for- 
mulation of transport in superlattices. We have found 
that the dynamics of the two-band model can show, in ad- 
dition to conventional Bloch oscillations, significant addi- 
tional structure: Zener resonances, stable plateaus, and 
band collapses. 
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FIG. 1. Energy spectrum for a two-band superlattice as 
a function of applied dc-field; continuous line corresponds 
to Eq.(0), while asterisks represent a model with no inter- 
band coupling {R — 0). The superlattice parameters are 
d = lOnm, A? = 0.8 x lO^^eV, = 0.92 x IQ-^eV, 
A'"' = Ag - Ag = 2.0 X 10~^eV, and R 
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FIG. 2. Time-dependence of p-{K = 0) for three different 
field values, top: eE = 0.9 x lO''; middle: 1.1 x 10^; and 
bottom: 1.02 x 10^ eV/m, respectively. For clarity, we have 
shifted the top and middle curve upwards by 0.4 and 0.2, 
respectively. The superlattice parameters are as in Fig. |l| 
and the unit on the time axis is 10'^?i/eV~ 4.14ps. 

FIG. 3. The negative of p™'" as a function of applied field. 
If — pll"" ~ 1, a strong band-to-band transfer is taking place 
('Zener resonance'). The dashed lines mark where the Wan- 
nier-Stark- ladder separation has local minima. 

FIG. 4. Fourier spectrum of p-{K = 0,t). A dot corre- 
sponds to each significant peak in the Fourier spectrum and 
the continuous lines are energy differences between the inter- 
penetrating Stark ladders. 

FIG. 5. Time-dependence of p_ {K = 0) for superlattice 
parameters AJ = = 1.8 x 10"^eV, A*"" = 2.0 x 10"^eV, 
and R = —0.9d. Units for the time-axis are as in Fig 2, and 
the field is _B = 0.232 x lO^V/m, corresponding to the eighth 
Zener resonance. 

FIG. 6. Quasienergy spectrum for the superlattice of Fig. 
1. We display the field-dependence of 17 fc-points, which are 
evenly distributed in the positive half of the Brillouin zone, 
k £ [0,n/d] (due to evenness of the quasienergy spectrum 
negative k's would be redundant). The vertical lines indicate 
where the band collapse would occur for non-interacting mini- 
bands. The modulation frequency is tJac = 1.5 x 10~^eV/ft, 
and the quasienergies are given in units of cjac. 



